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Abstract 

We present a possible explanation of the recently observed 511 keV 7-ray anomaly with a new 
"millicharged" fermion. The new fermion is light (©(MeV)) but has never been observed by any 
collider experiments mainly because of its tiny electromagnetic charge ee. We show that constraints 



,^ i from its relic density in the Universe and collider experiments allow a parameter range such that 

the 511 keV cosmic 7-ray emission from the galactic bulge may be due to positron production from 



this millicharged fermion. 
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The SPI/INTEGRAL observation of the very sharp 7-ray peak at 511 keV from the 
galactic bulge [l|, |2] needs an explanation of its origin. Most probably, it may come from 
the positronium decay. For this explanation of the positronium decay, a sufficient number 
of positrons are needed in the ffist place. The positron abundance in the galaxy can arise 
from several origins. 

Some obvious candidates are the astrophysical production mechanisms of positrons dis- 
cussed in |3| . However, these mechanisms through astrophysical sources such as black holes 
and supernovae turn out to be inappropriate to explain the intensity of the positron annihi- 
lation flux, especially in emission region, because the astrophysical sources like black holes 
and supernovae are expected to be more spread out than observed. Therefore, now most 
preferred interpretations of the 511 keV 7-rays rely on particle physics origins where new 
particles beyond the standard model (SM) are introduced.^ Usually a new particle in the 
mass range 1 — 100 MeV is introduced |5|.^ Let us call this new particle x- Recent anal- 
ysis including the internal bremsstrahlung radiation and in-flight annihilation gives more 
stringent mass bound for the light particle in MeV region: m < 3 — 4 MeV.^ The needed 
positron abundance may arise from the x decay and/or x~X annihilation to e^e~ . The new 
light particle should have negligible couplings to photon and Z boson; otherwise, it must 
have been observed at the LEP experiments. If the new particle is neutral under the gauge 
transformations of the SM as a heavy neutrino, it overdoses the Universe as noted by Lee 
and Weinberg 8|]. Thus, we exclude the neutrino possibility toward the origin of the 511 
keV line. This has led to a new particle, coupling to another gauge boson beyond the SM, 
e.g. as in Ref. [9|]. 

If another light ^(1) gauge boson, which will be called "exphoton, "^ beyond the SM exists. 



^ We have notieed a reeent elaim that the 511 keV hne distribution reported in the newest result from 
INTEGRAL seems to resemble the lopsided distribution of the "hard" low mass X-ray binaries (LMXBs) 
(low mass x-ray binaries with strong emission at E^ > 20 keV) [j|. However, quantitatively improved 
understanding of 511 keV gamma ray flux coming from LMXBs is required to see if the LMXB can fully 
account the anomaly. More observation would also be required for this issue. 

^ See also [6[ where O(IOO) GeV weakly interacting massive particles are considered. 

^ This constraint can be released by a factor of two by a possible ionization of the propagation medium [7[ . 

^ In the literature, the term "paraphoton" is commonly used. However, we use "exphoton" to emphasize 

the word "extra" which only directly couples to the "extra" matter field x ^'iid it is the gauge boson of 
the "extra" U{1). Moreover, this can show the fact that the extra Eg gauge group may contain exphoton 
in heterotic string models. 



most probably a kinetic mixing can exist via loop effects [10| between photon and exphoton 
without violating the charge conservation principle. After a proper diagonalization procedure 
of the kinetic energy terms, then the electromagnetic charge of x can be millicharged. In 
heterot. stnng models, the extra E', gauge group .nay coutaan the exphoton, leadurg to 
the kinetic mixing [11|. Indeed, an explicit model for this kind from string exists in the 
literature |l2l |. 

Very light ((9(eV)) millicharged particles with a sufficiently small charge are phenomeno- 
logically acceptable as studied in recent papers 13|]. On the other hand, the heavy mil- 
licharged particle idea as a dark matter (DM) candidate was suggested about 20 years 



ago 



14| and it has been revived recently 15|]. The intermediate 0(MeV) millicharged par 



tides has not been ruled out by observations in the previous study 16|] which, however, did 
not include the 511 keV line possibility. Earlier, the (9(MeV) millicharged particle effect 
on cosmic microwave background radiation was studied in the parameter region of expho- 



ton coupling constant (ocx = ^1^/^t^ ~ 0.1) [17|]. Here, we analyze the urgent problem of 



the (9(MeV) millicharged particles toward interpreting the 511 keV line within the limit 
provided by the DM constraint with reasonable exphoton coupling constants.^ 

Consider two Abelian gauge groups f/(l)QED and f/(l)cx-^ The kinetic mixing of f/(l)QED 
photon and ?7(l)cx exphoton is parameterized as 



C = -^F„,,F^^ 



flU-' 



iX.^x^'^ 
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where A^(X^) is the U{1)qed{U{1)cx) gauge boson and its field strength tensor is F^y{X^y). 
The kinetic mixing is parameterized by ^ which is generically allowed by the gauge invariance 
and the Lorentz symmetry. In the low-energy effective theory, ^ is considered to be a 
completely arbitrary parameter. An ultraviolet theory is expected to generate the kinetic 
mixing parameter ^ 10| . The usual diagonalization procedure of these kinetic terms leads 



^ The laboratory and cosmological bound of millicharged particles was studied sometime ago |l6| , but the 

study toward 511 keV line and the subcV mass range has not been included. 
^ One should note that the U{1) mixing in the observable and hidden sectors should be considered carefully. 

For the simple assumption of the charges given in Ref . 18| , x coupling with the full strength to the massive 



exphoton does not couple to the massless photon, or at least suppressed by e. Converting this argument, 
the massive Z-boson mixing with the massless exphoton gives the neutrino coupling to the exphoton 
suppressed by e. Thus, the very stringent supernova cooling constraint which gives a bound for the low- 
energy dark matter (m > 10 MeV) [19[ does not apply to our case since ux cross section is suppressed by 
e^ compared to that of |19|. 



to the relation, 







(2) 



and we obtain 



C = -\f,^F'^^ - \x,^X^\ (3) 

where the new field strengths are F^y and X^,^. Photon corresponds to A^ and exphoton 
corresponds to X^. If the exphoton is exactly massless, there exists an 50(2) symmetry 
in the A^ — X^ field space: A^ -^ cosOA^ + sin6'X^ and X^ -^ — sin 6*^4^ + cosOX^. Any 
physical observable, however, does not depend on 9. 

Using the above 5*0(2) symmetry, let us take the following simple interaction Lagrangian 
of a SM fermion, i.e. electron, with a photon in the original basis as 

C = i,{eQr)^K- (4) 

Note that in this basis there is no direct interaction between the electron and the hidden 
sector gauge boson X. If there exists a hidden sector Dirac fermion x with the ?7(l)cx charge 
Qx) its interaction with the hidden sector gauge boson is simply represented by 

/: = x(eexQx7^)x^^, (5) 

where Ccx can be different from e in general. In this case, there is also no direct interaction 
between the hidden fermion and the visible sector gauge boson A. We can recast the 
Lagrangian (jl]) in the transformed basis A and X, 

^ = ^(-7=^7^^^.. (6) 

Here, one notices that the standard model fermion has a coupling only to the visible sector 
gauge boson A even after changing the basis of the gauge bosons. However, the coupling 



constant e is modified to e/y 1 — ^^, and so the physical visible sector coupling e is defined 



as e = c/a/I — ^^. Similarly, we derive the following for x. 



C = xr\ eoxQx^M - ecx-7==gx^M 1 X- (7) 
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FIG. 1: The "millicharge" annihilation diagrams to, (a) e^e and (b-1) 27cx; (b-2) the 
brenisstrahlung diagram related to (b-1). The cross diagram in (b-1) is not shown. 

In this basis, the hidden sector matter field x iiow can couple to the visible sector gauge boson 



A with the coupling —e^^^/^/l — ^^. In terms of the aforementioned 50(2) symmetry, it 
simply means the mismatch between the gauge couplings of the electron and other fermions. 
Thus, we can set the physical hidden sector coupling e^x as eex = Ccx and we define the 
coupling of the field x to the visible sector gauge boson A, introducing the millicharge 
parameter e, as ee = —ecx^/^Jl — ^^. Note in general that e 7^ Ccx- Since ^ ^ ^e/cex is 
expected to be small, the condition ^ < 1 gives ac^/a > e^. From a fundamental theory, 
one can calculate the ratio ep:Je in principle, which is possible with the detail knowledge of 



the compactification radius 20[. Here, we simply take the ratio as a free parameter. 
For the cosmological study of Xi ^^ need the annihilation cross sections of DM: XX 



e e"" 



XX 



27ex, XX ~^ 77ex, and xX ~^ 11- The ratio for these cross sections is given by 
0-2^^^ : ae+e- ■ o'77ox ■ o-2'y — Oi^^ : E a : e aaex ■ £ a . (8) 



We noticed that the first two channels (depicted in Fig. [T]) are important and the last two 
channels are quite suppressed in the parameter region where e and ac^/a are quite small as 
is required by the observational data. If otc^ja > 0.01(0.1), the background diffuse gamma- 



ray fiux could be larger than 1(10)% of the 511 keV f 



ux, so the region is already excluded 



by the INTEGRAL and COMPTEL measurements [Tl, l2l| (See Fig. E]). As we will see 
below, XX ~^ 27ex channel ((b-1) in Fig. [T]) overwhelmingly dominates in the first two main 
channels of Fig. [H Then it seems that the gamma-ray fiux from the real bremsstrahlung 
((b-2) in Fig. [1]) could be of considerable amount. However, the bremsstrahlung cross section 



is suppressed by a factor of e'^a compared to that of diagram (b-1) in Fig. [H Thus, o"2^^ 

5 



brcm 



Q;(aex/«)^cre+e- and is negligible. The annihilation cross sections determine the relic density 
of the hidden sector fermion x- The process XX ~^ e~e'^ determines the flux of the eventual 
511 keV photons as well. Let us assume that the charge of the x particle is (0, Cex) in the 
basis of {A, X). The millicharge ee comes from the shift of the exphoton field in Eq. ([7]) and 
Cex is for the hidden sector f/(l)ex gauge interaction. 

The cross section for the process xX ~^ e~e~^, shown in Fig. [T](a), is given by 



XX^e e"' 
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(9) 



where Pi = \J\ — Arrif/s is the velocity of the particle-z and a = e^ /An. In the nonrelativistic 
regime, the approximation E ~ m-^ + i?7i^(fi.oi/2)^ makes sense and we obtain 

1 1 . m; 



(J 
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Now, the cosmologically interesting average of the cross section times velocity, {av)e-e+, 
becomes {av)e-e+ = fle-e+ +be-e+{'v'^) + C'((t'^)) where ae-e+ and 6e-e+ are given by Eq. ( TTOl) . 
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Similarly, for the process xX ~^ '^Icx shown in Fig. [T](b-1), we obtain 



da 



XX^27o: 



dcosO 






1 + 2plsm^9 - p^{2sm'^9 + cos 



p^cos^ey 



(12) 



where c^ex = egx/4vr. The total cross section is given by o";^^^27cx = /o '^i^'^^^)'^^- -'-^ ^^^^ 
case also, the cosmological average of the annihilation cross section times velocity, {o'v)2'y^^, 
is expressed in powers of v'^ as {av)2^^^ = a27ox + ^273x(^^) + ^ii'^'^)) where 027^^ = ^"cx/"^y 



and 627CX = ||o27cx- Again, we neglected the contributions from xX 
the smallness of e and Oc^/a. 

The relic density of a generic relic, X, can be expressed as 

2 1.07 X lO^GeV^^ xf 1 



77cx5 77 because of 



fiv/l 



8.77 X lO^^^MeV^^ 



5f^ (a + Sb/xp) 

Xp 1 

■^(a + 3b/ xp) 



(13) 
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where g^ is evaluated at the freeze-out temperature Tp, a and h are the velocity independent 
and dependent coefficients, respectively, andxp = mx/Tp — 17.2+ln((^/(^*)+ln(mx/GeV) + 



luy^aTp ~ 12 — 19 for particles in the MeV-GeV range (22[. We can approximate xp ~ 
11.6 + ln(?72x/MeV) for 1 MeV < nix ^ 100 MeV. Therefore, we can estimate the relic 
density of the millicharged particle, Xi &s 

^ ,9 ^■^ (11.6 + In m)m^ , , 

n^h^ ^ 1.60 X 10"^^ ^ ^ „ . 1/2 . TV' (14) 



(^)' + ^^(l-i) (l + 4 



where m = m^/MeV and we put g^ :^ 10.75 for 1 < Tp/MeV < 100.^ Finally, we can find a 
constraint for the mass m^ and the charge e of the millicharged DM and the hidden sector 



coupling Oex, based on the relic density of DM from the WMAP three- year results [23|. In 
Fig. [21 we present the excluded parameter space for typical DM masses (m^ = 1,3 and 10 
MeV) as the yellow shaded regions from our analysis of the DM relic density. 

In charged medium, photon can effectively obtain mass via the interaction with charged 
particles. Therefore, this effective mass should be smaller than the limit of photon mass. As 



a result, the Debye screening length in the DM medium around Earth Ad = y^x/^^^hc 
is required to be larger than the limit of the inverse photon mass J2J, |25|. Putting n^ = 
Pxl'^x — 0.3GeV/cm'^ x ^^/(fiDM'^x) ^"^^ ^dm — 0.23, we finally get the simple relation 
^ > 2825. One should note that the relic density of x is essentially proportional to m? so 
that the Debye screening length is not sensitive to the mass. The lower right corner from the 
central region (colored by pink) is excluded by this constraint. Interestingly, m^ > 3 MeV 
does not have the parameter space which can fully accommodate the dark matter density 
^DM ^ 0.23. 

The line e = oic^ja corresponds to the line of equal couplings that divides where the 
diagrams (a) and (b-1) in Fig. [T] dominate: in the upper part of the line the process XX ~^ '^Ic^ 
and in the lower part the process xX ~^ e~e~^ dominate. In addition, we show the allowed 
range of e for typical DM masses {m^ = 1,3 and 10 MeV) as the (gray shaded) vertical bands, 
which will be obtained from the following analysis of the 511 keV 7-ray flux constraint. For 
example, if m^ = 3 MeV, the middle (gray shaded) vertical band for e in the upper left 
corner is allowed. The smallness requirement of ^ is buried in the Debye screening length 



constraint. The study of 17| is buried in the lower right corner around ^ = 1. As can 



^ In this step, we use the total annihilation cross section, i.e. a = ag-g+ + a27„x ^^'^ ^ — ^e-e+ + ^: 



27e: 
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FIG. 2: The plot for a^^/a versus e. The lower left corner (yellow shaded regions) is excluded 
by the DM relic density constraint: the lines correspond to Q^h? = 0.11 and m^ = 1,3, and 
10 MeV, respectively. The vertical bands (gray shaded regions) are the allowed range of e that 
will be given by the 511 keV 7-ray flux constraint analysis: the regions correspond to m-^ = 1,3, 
and 10 MeV, respectively. The region excluded by the Debye screening is shown from the central 
(pink shaded) region to the lower right corner marked by Debye screening. The (green) region 
Oex/a > 0.01(0.1) is excluded since more than 1(10)% diffuse gamma ray flux compared to the 511 
keV flux is expected. 



be seen from the figure, a significant region is excluded. However, we note that there still 
remains an available space. 

The observed flux of dark matter annihilation products can be obtained by integrating 
the density squared along the line of sight as 

dN, 1 



$i(?/^,E) = av- 



dsp {r{s,^jj)) 



(15) 



dE 4:7Tml^ Jiine of sight 

where p{r) is the mass density of the DM, a is the DM annihilation cross section, v is the 

velocity, dNi/dE is the spectrum of secondary particles of species i, and s is the coordinate 



running along the line of sight, in a direction making an angle, ib^ 



rom the direction of the 



galactic center. It is convenient to introduce the quantity J{'ip) 261]: 

jfw,) = - I Z. ) / dso^iris^)) 

8.5 kpc V0.3GeV/cmV Ane of sight ' "^ "- ' "^ " 



(16) 



by which the expression in eq. (fT5|) can be separated to "halo profile depending" factors and 
"particle physics depending" factors as 

dNi fav\ /lMeV\% 



, „ ,_.., _ (17) 

dE \pb/ \ mDM / 

where J(Af2) is defined as the average of J^^j) over a spherical region of sohd angle, AO, 



<l>i(Afi, ^) ~ 5.6 — '-[ — ] — J (Afi) Afi cm-^s-i 

dE Vpb/ V "^DM / 



centered on t/^ = 



22|. 



If the mass of the DM particle is less than the muon mass, the low velocity annihila- 
tions can produce electron-positron pairs. Most positrons lose energy through their interac- 
tions with the inter stellar medium (ISM) and bremsstrahlung radiation and go rest. Then 

n 

positron annihilation takes place via the positronium formation (~ 96.7 ± 2.2%) [2] and 
partly via the direct annihilation into two 511 keV gamma-rays. Only 25% of the time, a 
singlet positronium state decaying to two 511 keV photons is formed while 75% of the time, 
a triplet state decaying to three continuum photons is formed. This means that the 511 
keV photon emission occurs only by a quarter of the total positron production through DM 
annihilation. After taking all this into account, the flux of 511 keV 7-rays from the galactic 
center can be given as 

{(yv\ /iMeVV 



$7,511 ^ 



0,275 X 5,6 ^-j 



m. 



■; 



-2„-l 



J(Afi)Aficm-^s 



(18) 



where Afi is the observed solid angle toward the direction of the galactic center. 

The observed 7-ray profile has a full width at half maximum of ~ 6° with a 4° — 9° 2a 
confidence interval and the flux $7,511 ^ (1.02 ± 0.10) x 10"^ ph cm~^ s~^ |ll;l2i- Thus, we 
consider a solid angle of 0.0086 sr, corresponding to a 6° diameter circle. In this model, 
positron is produced from the process XX ~^ e~e~^. Therefore, we can find the charge e of 
the millicharged DM as a function of its mass m^ from the resultant cross section (crt')e-e+ 
for this process and Eq. flTSl) . The relation is given by 



£ ~ 1.0 X 10"^ 



2 

m 



Vj 



m 
m' 



-1/4 



1 + 



m2_ 
2ml 



-1/2 



(19) 



where m = m^/MeV. To estimate the required parameter space, we use the width of the 
observed distribution J(0.0086 sr) ~ 50 — 500, approximately corresponding to 7 ~ 0.6 — 1.2 
essentially following the approach of Ref. 27|.^ 



If the main source of 511 keV 7-rays from galactic bulge is from the DM annihilation, the observed 
distribution of 511 keV emission line would constrain the shape of the DM halo profile because DM 
annihilation rate is proportional to the DM density squared. 
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There already exist various bounds from experimental and observational results, which are 



summarized in 16|] . Among them, the limit from the millicharged particle search experiment 



at SLAG 



is relevant to the mass-charge parameter space, which is considered in this 



analysis. In principle, the DM can contribute to the anomalous magnetic moment 29|] , but it 
can only occur at the two-loop level with an additional e^ suppression factor. The expected 
recoil energy by the DM-nucleon scattering is too small to be measured by the existing or 
near-future experiments because of the lightness of the proposed DM candidate. 

The result from the study of the 511 keV 7-ray flux and the SLAG experiment is presented 
in Fig. [3] in the e — m^ space. Even after taking into account the SLAG bound for the 
millicharged particle, a large parameter region is still remaining. Recent analysis such as 
the internal bremsstrahlung radiation and in-flight annihilation gives strong mass bound for 
the light dark matter in MeV region: ?Ti < 3 — 4 MeV.^ Therefore, the lower left corner is 
magnified. In the allowed parameter region (1 ^ «cx/ct > s\ the relic density of DM is 
essentially determined by XX ~* "^Icx- However, the observed 511 keV photon flux is mostly 
explained by xX ~^ e~e+. In this respect, the difficulty of explaining both quantities in 
Ref. J27| is easily avoided in our model. 

One final comment is about the spontaneously broken t/(l)cx symmetry which results 
in the nonvanishing exphoton mass. In this case, the electrically charged particles such as 
electron and proton can couple to the exphoton though the hidden fermion (x) does not 



directly couple to the on-shell photon 18|]. In principle, this case can be also relevant to our 
DM problem and the related 511 keV photon line. Theoretically, spontaneous symmetry 
breaking generally gives finite ranges of parameter space both for massless and massive 
exphotons and hence our study on massless exphoton covers a finite range of the parameter 
space. In the future, we would like to discuss the cosmology of (9(MeV) exphoton. 

In conclusion, we presented an allowed parameter range of a new millicharged particle 
X with (9(MeV) mass toward a possible solution to the recently observed 511 keV cosmic 
7-ray anomaly. It couples to photon with a "milli" electric charge strength, ee. In the mass 
range of m^ < 1 MeV, the millicharged particle can constitute a sizable (> 10%) portion 
of the DM content of the Universe but might have escaped detection so far in any collider 



^ As already stated in the beginning, this constraint can be reduced by a factor of two by a possible 
ionization of the medium 
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FIG. 3: The plot for e versus m^. The dark (blue) shaded region is excluded by the SLAG search 
of millicharged particles 



28l |. We plot the Navarro- Frenk- White profile [30|] model line (dotted) 
using the fitting parameter {R = 20 kpc, po = 0.347 GeV/cm'^) and the lightly shaded (yellow) 
region for the uncertainty range J(0.0086 sr) ~ 50 — 500. After considering the recently given 
strong constraint on the light dark matter mass [71], the allowed region is m^^ < 4 MeV (bottom). 

experiments basically because of its tiny electric charge. This millicharged particle may 
arise in a more fundamental theory such as string as an interplay between the observable 
and hidden sectors. 
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